An Improved Nonperturbative Method for Studying Two-dimensional Vortex Liquids 
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We show how a systematic improvement can be made on the nonperturbative parquet approxima- 
tion method which was previously used to study the effect of thermal fluctuations in vortex liquids in 
high-temperature superconductors. This is achieved by including an infinite subset of Feynman di- 
agrams contributing to the renormalized four-point vertex function of the Ginzburg-Landau model, 
which was omitted in the original approximation. We find that the growing crystalline order in 
the vortex liquid is more pronounced in the improved approximation. In particular, the second and 
third peaks in the liquid structure factor, which appeared as one peak in the original approximation, 
are now resolved. 

PACS numbers: 74.20.De, 74.60.-w 
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I. INTRODUCTION 

Thermal fluctuations play a more important role in 
high-temperature superconductors than in conventional 
materials because of strong anisotropy, high tempera- 
ture, and short coherence lengtbA Therefore, a high- 
temperature superconductor in a magnetic field is be- 
lieved to be in a vortex liquid phase resulting from the 
melting of the Abrikosov vortex lattice predicted by the 
mean field theory^. For bulk materials, experiments^^ 
detect well below the upper critical field H C 2 sharp drops 
in resistivity and steps in the magnetization and the spe- 
cific heat, which are interpreted when the strength of 
disorder is weak as being due to a vortex liquid under- 
going a first-order phase transition into presumably the 
Abrikosov lattice as the temperature is lowered. How- 
ever, the situation is more complicated because the first- 
order transition disappears at both high and low mag- 
netic fieldo 3 i 6 i 7 i 8 . 

For two-dimensional superconducting films, the ques- 
tion whether the vortex liquid undergoes a phase tran- 
sition at all into the low-temperature ordered state is 
still controversial. The perturbation expansions around 
the high-temperature liquid state using the Ginzburg- 
Landau (GL) model within the lowest Landau level 
(LLL) approximation have been carried out for evi- 
dence of the transition to the ordered phase at low 
temperatures&iii. A weak first-order transition is de- 
tected in two-dimensional vortex liquids from numerical 
simulations 1 -. However, numerical simulation results de- 
pend strongly on the boundary conditions as the Montc- 
Carlo simulation performed on a spherical geometry 12 
shows no sign of a finite-temperature transition. 

As a nonperturbative analytic (i. e. nonsimulational) 
approach, the parquet approximation has been success- 
fully applied to a two-dimensional vortex systemic and 
also to vortices in a layered superconductor^^, ft is 
free from any finite-size effect perpendicular to the field 
direction and sophisticated enough to capture the grow- 
ing crystalline order developing in the vortex liquid as 
the temperature is lowered. The parquet approximation 
deals with the renormalized four-point function of the 



vortex system which is obtained by summing an infinite 
subset of Feynman diagrams, the so-called parquet di- 
agrams. Although the parquet diagrams seem to form 
a minimal set of Feynman diagrams which capture the 
growing crystalline order properly, there is no a priori 
reason to neglect the non-parquet diagrams. It is be- 
cause there is no apparent small parameter associated 
with the non-parquet diagrams. In this paper, we present 
a first attempt to go beyond the parquet approximation. 
We show how one can include systematically the non- 
parquet diagrams into the existing nonperturbative cal- 
culations. We do this by devising a way to take into 
account yet another infinite subset of Feynman diagrams 
which arc omitted in the previous calculations. This pro- 
cedure is similar to what one does in the integral equation 
approach to ordinary liquids^ where integral equations 
such as the hypernetted chain equation are improved by 
adding an appropriate set of diagrams. We shall see that 
the growing crystalline order developing in the vortex 
liquid as the temperature decreases is more pronounced 
when the non-parquet diagrams are included, although 
no finite temperature phase transition is detected as in 
the previous studies^. 

In the next section, we briefly set up the parquet dia- 
gram decomposition method for the two-dimensional vor- 
tex liquid systems. In Sec. IIIII we show how the non- 
parquet contributions can be included in the formulation. 
The improved integral equations which include these new 
diagrams are solved numerically in Sec. II VI We conclude 
with discussion in Sec. 



II. PARQUET GRAPH RESUMMATION 
METHOD 



The parquet graph resummation method for vortex 
liquids^ is based on the LLL approximation of the GL 
model for a superconductor in a perpendicular magnetic 
field. For a superconducting film, the GL free energy 
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with the order parameter denoted by is given by 
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where a, /3, and m are phenomenological parameters and 
e* = 2e. We take B = V x A as constant and uniform, 
and use the LLL approximation, which is believed to be 
valid over a large portion of the vortex-liquid region— In 
the symmetric gauge, where A = ^(— y,x, 0), the LLL 
wavefunction is given by ^]/ LLL (r) = exp(— fj, 2 \z\ 2 /i)<f>(z) 
where /i 2 = e*B/hc and <fi(z) is an arbitrary analytic 
function of z = x + iy. In the LLL approximation, the 
free energy becomes 
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where an = oi + e* Bh/2cm changes sign crossing the 
upper critical field line H C 2(T). The effect of ther- 
mal fluctuations in the two-dimensional vortex liquid 
systems is determined by the partition function Z = 
J T> (pT) <fi* ex-p(—F[(f>, </>*]). In this section, we set up the 
parquet diagram resummation method and calculate var- 
ious correlation functions with respect to the partition 
function. One can develop the standard perturbation 
theory from the given partition function. The bare prop- 
agators are given by 
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Since the magnetic length is the only length scale 
perpendicular to the field direction which appears in the 
propagator—, the fully renormalized propagator can also 
be written in the same way as Eq. (J3J with the renor- 
malized a a replacing the bare ajj- The renormalized an 
is determined self-consistently in the parquet approxima- 
tion as will be seen later. 

The main quantity one calculates in the parquet ap- 
proximation is the renormalized connected four-point 
function ((f)* (z*)<f>* (z*)(j)(z^)(j){zi)) c . An important fea- 
ture of the LLL approximation is that this renormalized 
quantity can be completely described by a single vertex 
function r(k)a*, where the momentum k corresponds to 
the two-dimensional space perpendicular to the magnetic 
field. In general this quantity can be written as 
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FIG. 1: A graphical representation of the Bethe-Salpeter 
equations for the reducible parts IL. (See Eq. J7J.) 



Note that to the lowest order, T(k) = r.g(k) = 
exp(— \k\ 2 /2fj, 2 ) is the bare vertex. 

The parquet approximation has been widely used— in 
many branches of many-body physics. For the vortex liq- 
uids, we make a resummation over all parquet diagrams 
by first noting that the contributions to T can be decom- 
posed into the totally irreducible part denoted by R and 
the reducible part. The reducible part in turn can be 
written as the sum of three parts IL, (i = 1, 2, 3) repre- 
senting the contributions from three different channels as 
shown in Fig. ^ (A detailed discussion of the diagram- 
matic decomposition can be found in Ref. El) We have 



r(k) = J R(k) + ^n i (k). 



(5) 



Each reducible vertex IL is composed of an irreducible 
vertex A; where 



A J (k) = J R(k) + ^n,(k) 



(6) 



and the renormalized T via the following Bethe-Salpeter 
equations (see Fig.^l: 



n 1 (k) = -x[A 1 or] (k), 
n 2 (k) = -2sA 2 (k)r(k), 

n 3 (k) = -2:r [A 3 *T] (k), 



(7a) 
(7b) 
(7c) 



where the operation o between two arbitrary functions 
/(k) and g(k) is defined by 
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and / * g is just the convolution without the cosine fac- 
tor. The parquet approximation corresponds to keeping 
only the bare vertex function I\s(k) in the totally ir- 
reducible vertex i?(k), and neglecting all the higher or- 
der non-parquet diagrams. The lowest order of the di- 
agrams neglected in this approximation is 0(/3 4 ). The 
main point of this paper, which will be discussed in the 
next section, is to find a way to incorporate systemati- 
cally the diagrams neglected in the parquet approxima- 
tion. In Eq. 0, we have used the dimensionless pa- 
rameter x = /i 2 f3 /2ita 2 R , which originated from the two 
renormalized propagator lines in IL,. This parameter x 
is determined self-consistently from the following exact 
Dyson equation, which relates the renormalized propa- 
gator and the renormalized quartic vertex: 
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where ar = c%h -J^tt/ f3fi 2 is the dimensionless tem- 
perature. We can show that ax is proportional to 
-(1 -t- hj/ith) 1 / 2 , where t = T/T cQ and h = H/H c2 (0) 
with T c0 and H C 2(0) being the critical temperature at 
zero field and the upper critical field at zero tempera- 
ture, respectively. Note that Eqs. ©, ©, and © 
form a closed set of equations for T(k) for given ar and 
R(k). Note also that this set of equations are exact re- 
lations for the vertex function T(k). Recall that in the 
parquet approximation i?(k) is approximated as the bare 
vertex r#(k) . 

Using the solutions to the above equations we can 
calculate several interesting physical quantities. Among 
them, we focus on the structure factor, which is a mea- 
sure of the correlation between vortices in a vortex liquid. 
It is calculated from 

X(r-r') = (|*(r)| 2 |*(r')| 2 > - (|*(r)| 2 ) (*(r')| 2 ). (9) 

The structure factor A(k) used in this paper is then de- 
fined as 



A(k) = ( 



y^Pj e k2 / 2 ^ J d 2 Re lk R X (R). (10) 



By joining two external legs of the four-point correlation 
functions in Q), we obtain a simple relation, 



A(k) = 1 - 2xT(k). 



(11) 



The above coupled integral equations for T(k) can be 
solved numerically. In Ref. I13L the parquet equa- 
tions were solved for two-dimensional vortex liquids with 
and without quenched impurities. The parquet equa- 
tions can also be solved for the vortices in layered 
superconductors^. 




FIG. 2: The leading-order non-parquet diagrams. Note that 
a general vertex F p (k) is used on each of four vertices. There 
are two more diagrams obtained by exchanging 23 and 24. 



III. NON-PARQUET CONTRIBUTIONS 

In this section we discuss how the contributions from 
the non-parquet diagrams can be included into the corre- 
lation functions of the two-dimensional vortex liquid. Re- 
call that, in the parquet approximation, all the diagrams 
contributing to the totally irreducible vertex i?(k) are ne- 
glected except the bare vertex diagram. The lowest order 
of the non-parquet diagrams is 0(/3 4 ) as shown in Fig. [21 
A straightforward way to proceed would be to calculate 
the non-parquet diagrams term by term starting from 
the fourth-order diagram. This procedure will generate 
a perturbation series for i?(k). To extract nonperturba- 
tive information for i?(k) and for the structure factor, 
one must then apply a resummation method such as the 
Pade approximation to the perturbation series. Such nu- 
merical resummation procedures, however, are knownifi 
to be less effective in capturing the growing crystalline 
order in the vortex liquid than the integral equation ap- 
proach like the parquet approximation. 

In this paper, we take a different route by incorporat- 
ing an infinite subset of the non-parquet diagrams sys- 
tematically into the integral equations discussed in the 
previous section. We first evaluate the diagrams repre- 
sented by the lowest order non-parquet diagram shown 
in Fig. |2 Note that, in place of the bare vertices, we 
use a general vertex T p (k) to be specified later. For any 
vertex r p (k), these diagrams certainly represent a subset 
of Feynman diagrams that is not considered in the par- 
quet approximation i. e. that cannot be decomposed as 
in Fig. ^ We find that there are two distinct diagrams 
as shown in Fig. [2 which turn out to give the same con- 
tribution denoted here by J(k). By explicitly evaluating 
these diagrams using a general vertex r p (k), we obtain 
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There are also contributions from the diagrams obtained 
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FIG. 3: A higher-order non-parquet diagram which can be 
included using the present method 

by exchanging 2:3 and 24 in Fig. which are related to 
J(k) via the hat-operation defined by 

? 2vr f d 2 p fk x p y -k y p x \ 

J(k) = ujJwr J(p) cos I — 7 — ) • (13) 

To include the contributions from the new set of dia- 
grams, we take in Eqs. JSJ and (JSJ 

i?(k) = r B (k) + j(k) + j(k). (w) 

Now we must specify what r p (k) is in Eq. I|12|) to close 
the self-consistent equations. In this paper, we take r p (k) 
as the sum of all the parquet diagrams, that is the so- 
lution of the parquet equations ©, © and iQ when 
R(k) is just r#(k) for given parameter x. Then the non- 
parquet diagrams in Fig. [21 represent an infinite subset 
of Feynman diagrams for which the skeleton diagram in 
Fig.0contain all the parquet diagrams on each of the four 
vertices. The inclusion of these contributions is achieved 
by solving again Eqs. ©, © and but with Eq. (JTU). 
There are of course other choices for T p than the present 
one, which will be discussed in Sec. [V] We find that 
the present scheme is relatively easy to implement nu- 
merically and to generalize to the next order. To get an 
improvement of the present approximation, we just need 
to calculate the contribution from the next order skele- 
ton diagram shown in Fig. [3| and include it into Eq. I|14|) 
in a similar way. 

IV. RESULTS 

As explained in the previous section, we solve numer- 
ically the self-consistent equations I© and Q to- 
gether with ifHI) for T(k). The vertex function r p (k) 
used in (|14|l is obtained by solving once Eqs. ©, © and 
(0) with i?(k) = Ts(k). All these numerical procedures 
are performed for given value of x. The temperature pa- 
rameter ax is determined from Eq. (jSJ). Since we are 
considering the vortex liquid phase, we only consider a 
rotationally symmetric case where all the functions de- 
pend on a dimensionless momentum K with K = k//z. 



The integral equations are solved numerically by itera- 
tion starting from an appropriate choice of initial T(k). 
The convergence of the iteration can be improved when 
the solution at slightly smaller value of x is used as the 
initial choice. We obtain the solutions for given values 
of x up to x = 100, which is the largest value we con- 
sidered. It corresponds approximately to ay = —11.9 
(see Fig. 0}. As we go down to lower temperatures, we 
have to increase the fc-space cutoff to accommodate the 
peaks appearing at large k in the structure factor and de- 
crease the grid size at the same time to capture the sharp 
first peak (see Figs. 03 and©. The number of iteration 
needed to get a convergence increases as the temperature 
is lowered. At the lowest temperature we considered, we 
needed about 1000 iterations to get a convergence. All 
these factors limit the temperature range where the nu- 
merical solution can be obtained. 

Figure 01 shows the relationship between the renormal- 
ized propagator ^px ~ a^ 1 and the temperature parame- 
ter ctT- The values are compared with those obtained 
when only the parquet diagrams are considered. We 
can see that there are essentially no difference in the 
two cases. Note that in the present analysis, we have 
extended the previous parquet approximation results^ 
to lower temperatures olt — —11.9. The dimensionless 
renormalized propagator *Jx is directly proportional to 
the thermodynamic quantities like the magnetization and 
the entropy of the vortex liquid system. Another interest- 
ing thermodynamic quantity is the generalized Abrikosov 
ratio defined by /3 A (x) = (|#(r)| 4 )/[(|W(r)| 2 )] 2 , where 
the bar denotes the spatial average. As shown in Ref. IT^l 
this quantity decreases as the temperature is lowered, 
and approaches 1.16, the value for a triangular lattice 
as x — ► 00 or ctT — * —00. It is related to T(k) in 
such a way that the relation © can be rewritten as 
aT = (1 — x/3a(x))/ ' \fx. Therefore, the generalized 
Abrikosov ratio even in the presence of the non-parquet 
contributions shows the similar behavior to the result 
obtained in the parquet approximation^. From these 
results, we may conclude that non-parquet contributions 
make little difference on the thermodynamic quantities. 
As we will see below, however, effects of non-parquet di- 
agrams appear in the crystalline order developing in the 
vortex liquid. 

We calculate the structure factors 111|) for the two- 
dimensional vortex liquid at various temperatures as 
shown in Figs. [5] and They are compared with those 
obtained in the parquet approximation. The structure 
factors for both cases look almost the same down to low 
temperatures ar ~ — 7. As the temperature is lowered 
further, however, we find that the first peak becomes 
slightly larger and sharper when the non-parquet dia- 
grams are included (see Fig. ©. This trend continues 
further down to lower temperatures. When the tem- 
perature is lowered below ax — —10, the second peak 
begins to split into two peaks (see Fig. ©. This can 
be interpreted as the non-parquet contributions captur- 
ing the growing crystalline order in the vortex liquids 
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FIG. 4: The renormalized propagator \fx ~ 1/ajj as a func- 
tion of the temperature parameter qt. The filled triangles 
and the circles are obtained with and without the non-parquet 
diagrams, respectively. The dashed line is a guide for the eye. 



more effectively. The peaks developing in the struc- 
ture factor correspond to the positions of the recipro- 
cal lattice vectors (RLV) of the triangular lattice. The 
RLV G can be represented in a dimensionless form as 
G//i = Go (mrj,n — mQ using a set of integers m and 
n, where, for the triangular lattice, 77 = \/3/2, £ = 1/2 
and Go = \j2n jr\ ~ 2.694. Therefore the lengths of the 
RLV can be grouped into |G|//x = c^Go, % = 1,2,3. . ., 
where c\ = 1, c% = v3j C3 = 2, C4 = v7> C5 = 3, etc. The 
first peak in the structure factor is located near K = Gq. 
Comparing the results in Fig. HJ1 for the cases with and 
without the non-parquet diagrams, we can see that when 
the non-parquet contributions are included, the first peak 
is closer to its expected position. Since, for a triangular 
lattice, the second and third sets of RLV are relatively 
closely spaced, they appear as one peak within the par- 
quet approximation. Figure shows that the resolved 
peaks are located near the expected positions of the RLV 
when the non-parquet contributions are included. The 
situation is similar for the closely-spaced fourth and fifth 
peaks. Although the peaks are not resolved at the mini- 
mum temperature we have studied, we can clearly see a 
tendency compared to the parquet results. 



V. DISCUSSION AND CONCLUSION 

In summary we have generalized the parquet approxi- 
mation for the two-dimensional vortex liquid systems by 
demonstrating that the non-parquet contributions can be 
included systematically into the nonperturbative calcu- 
lation of the correlation functions. The crystalline order 
developing in the vortex liquid as the temperature is low- 
ered is captured more effectively with the inclusion of the 
non-parquet diagrams in the sense that the peaks in the 



FIG. 5: The structure factor A of the two-dimensional vortex 
liquid as a function of the dimensionless momentum K — 
at qt — —7.66. The solid line is from the calculation 
including the non-parquet diagrams, while the dashed lines is 
that of the parquet approximation. The arrows indicate the 
positions of the RLV of the triangular lattice. 
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FIG. 6: The structure factor at qt — —11.9. The solid and 
dashed lines are the same as in Fig. The expected position 
of the first peak from the RLV is presented as a dotted line 
for a guide for the eye. 

structure factor corresponding to the RLV of the trian- 
gular lattice become resolved. 

In Ref. ll3l within the parquet approximation, no finite 
temperature phase transition has been observed in the 
two-dimensional vortex liquid as the temperature is low- 
ered. The length scale characterizing the growing crys- 
talline order in the vortex liquid was determined from 
the width of the first peak in the structure factor. It 
was shown that this length scale grows like y/x ~ \olt\ in 
the low temperature limit. Since the first peaks in the 
structure factor change only slightly with the inclusion of 
the non-parquet diagrams, these conclusions drawn from 
the parquet approximation remain valid. The detailed 
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structure factor has a more accurate form in the present 
generalization. 

As mentioned in Sec. IIIII there are other possible 
generalizations of the parquet approximation than the 
present one. They correspond to different choices for the 
vertex r p (k) used in Eq. El We can, for example, take 
r p equal to the full vertex T(k) not just the sum of all 
the parquet diagrams. Obviously the diagrams in Fig. [21 
in this scheme contain more diagrams than in the present 
one. We have attempted to solve numerically the set of 
equations ©, ©, and JT3J with (JT^I when T p = T. 
At relatively high temperatures, we find there is very lit- 
tle difference in correlations of the vortex liquid between 
this scheme and the present one. But the numerical it- 



eration in this case involves repeated evaluation of the 
integral in Eq. 1121 at every step of the iteration, which 
considerably slows down the whole calculation. We find 
that it is not practical to use this approximation below 
ax — —3. We believe that the effect of the non-parquet 
diagrams is already captured in the present approxima- 
tion scheme. 
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